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A NEW FORMULATION OF THE CONSERVATION 
EQUATIONS OF FLUID DYNAMICS 


By M. Vinokur 

Ames Research Center, NASA 
Moffett Field, Calif, 94035, U.S.A. 
and 

Department of Mechanical Engineering, University of Santa Clara 
Santa Clara, Calif. 95053, U.S.A. 


SUMMARY 

The computation of time-dependent flows has inspired a new, higher- 
dimensional formulation of the conservation equations of fluid dynamics in 
which time is treated as a fourth coordinate. The formulation is derived 
for a constant-density flow, and then extended to a variable-density flow 
by introducing a fifth, fictitious coordinate. This new coordinate can 
also act as a source coordinate, so that external source terms can be 
included. The analysis is carried out for both incompressible, stratified 
flow, and compressible equilibrium flow. The results are then extended to 
non-equilibrium and magnetohydrodynamic flows. Several applications of the 
new formulation to the computation of time-dependent flows are discussed. 
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1 . INTRODUCTION 

An important concept in computational fluid dynamics is the conservation- 
law form of the flow equations [1], These are differential forms of the basic 
integral relations expressing the conservation of mass, momentum, and energy. 
Since the latter are valid in the presence of discontinuities, appropriate 
methods of differencing the differential equations will result in the correct 
jump conditions being satisfied for the computed, smeared-out discontinuities. 
In Cartesian coordinates, the equations consist only of sums of derivatives - 
called a strong conservation- law form [2]. One can then difference the equa- 
tions so as to assure exact numerical conservation of mass, momentum, and 
energy for the total flow region. 

If one employs curvilinear coordinates, in order to satisfy conditions 
on general boundaries, undifferentiated terms due to the derivatives of the 
base vectors appear in the scalar decomposition of the momentum equation. 

This weak cons ervation= law form [2] prevents the achievement of overall 
numerical conservation of mass, momentum, and energy, and can thus lead to 
loss of accuracy in "capturing” sharply curved discontinuities. To avoid 
this, Anderson, Preiser, and Rubin [3] introduced Cartesian base vectors and 
achieved a scalar decomposition in the strong conservation-law form. A 
similar idea was used by MacCormack and Paul lay [4], although they employed 
the integral formulation of the equations. Calculations using this latter 
formulation have been carried out by Rizzi and Inouye [5] for blunt-body 
flows, and by Rizzi, Klavins, and MacCormack [6] for supersonic flows. 

Recently, the author [2] proposed a new, simpler, local scalar decomposition 
in the strong conservation- law form, using the curvilinear base vectors at 
the fixed, computational mesh points. 

There are two important classes of time -dependent problems in which it 
is useful to introduce new independent variables which depend on both time 
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and the spatial coordinates. One class involves time-varying boundaries such 
as flexible walls! detached shock waves, jet boundaries, etc. If the un- 
steady motion of these boundaries is coupled strongly to the time-dependent 
phenomena being studied, coordinate transformations making these boundaries 
coordinate surfaces will greatly simplify the satisfying of boundary condi- 
tions. The other class of problems is flow past (or inside) a rigid body 
undergoing arbitrary translational or rotational motion. It is then useful 
to employ a non-inertial reference frame fixed with respect to the body. 

In both cases, the time- dependent coordinate transformations result in 
undifferentiated "source" terms and a weak conservation law. The situation 
here is analogous to that discussed in the previous paragraph, where the 
replacement of Cartesian by curvilinear coordinates led to a weak conserva- 
tion law. It is therefore natural to seek a strong conservation- law form of 
the equations resulting from the time-dependent coordinate transformations, 
analogous to the strong conservation- law form derived by the author for 
curvilinear coordinates. If one does, one quickly finds that one must 
introduce a base vector associated with the time coordinate, and that the 
continuity and energy equations must be coupled with the momentum equation 
in a single conservation law. The analogy would then also require a fifth, 
fictitious coordinate (and corresponding base vector). The time-dependent 
coordinate transformation then becomes a simple coordinate transformation 
in the higher-dimensional space, and the five equations in strong conserva- 
tion form are given by the local scalar decomposition of a single vector 
equation. 

The author has developed a higher-dimensional formulation of the 
conservation equations of fluid dynamics, motivated by the above considera- 
tions, and has applied it in Ref. [2] to extend the strong-conservation form 
of the equations to time- dependent coordinate transformations. The purpose 
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of the present paper is twofold. One is to give a detailed derivation of the 
formulation (not presented in Ref. [2]), and where possible, to give physical 
justification for the mathematical assumptions that are made. The other 
purpose is to generalize the results to a wider class of flow equations, in 
the hope that the new formulation will find other fields of application. 

While the primary application is to generalized coordinates in compres- 
sible flow, in order to present the ideas as simply as possible the formulation 
is first derived in ^2 for orthogonal curvilinear coordinates in a 

constant-density flow. The flow equations are shown to result from a single 
four-dimensional equation in Euclidean space-time. The components in the 
spatial direction yield the momentum equation, while the component in the 
time direction gives the continuity equation. (By contrast, in the space- 
time of special relativity, time components are related to an energy equation). 
A physical interpretation of the time component of velocity is proposed. 
Finally, an examination of steady-state flow shows that the continuity equation 
can still be considered as the component of a four-dimensional equation in a 
fourth direction, even though the fourth coordinate is ignorable. 

The ideas in the previous sentence are used to consider variable-density 
flow in ^3. An ignorable, fictitious fifth coordinate is introduced, 

whose directions yield the additional conservation equation which is required. 
This coordinate can also be employed to include real source terms in the 
analysis. Both compressible equilibrium flow and incompressible stratified 
flow are treated. In ^4 the analysis is extended to non- equilibrium 

and magnetohydro dynamic flows. The restriction to orthogonal coordinates 
is also removed in this section, where results are given in general, curvi- 
linear coordinates. The discussion in ^5 centers on thp closure 

problem of completing the set of conservation equations by auxiliary 
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equations. As an illustration of the use of the formulation presented in 
this paper, this section concludes with a derivation of the scalar decomposi- 
tion in strong conservation- law form of the ideal gas equations for a 
general class of time-dependent coordinate transformations. A brief 
discussion of the conservation laws in a non-inertial reference frame is 
also presented. 


2. CONSERVATION EQUATION FOR CONSTANT -DENSITY FLOW 
2 . 1 Three-Dimensional Formulation 

Consider the non-relativistic motion of a fluid with constant density 
p. The continuity equation 

div (p2f) = 0 (2.1) 

expresses the conservation of mass, where u is the velocity vector. The 
momentum equation relates the time rate of change of momentum per unit 
volume to the net convection of momentum and the net force acting on the 
fluid contained in that volume. The surface force can be expressed as the 
divergence of a stress tensor. If the pressure p characterizes the state 
of stress existing in a fluid in uniform motion, the stress tensor for a 
fluid in non-uniform motion can be written as 
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S - pi , 

where the symmetric tensor ? is called the extra stress tensor or viscous 
stress tensor, and t is the identity tensor. In an inertial reference 
frame, the external body force in the absence of electromagnetic effects 
is due to the gravity. (An external force due to a magnetic field is dis- 


cussed in ^4.2.) Introducing the gravitational potential <f>, the body force 
per unit mass is 


f = - grad <f> = - div (01) 

• 

(2.2) 

If one defines a generalized stress tensor 



P = S- (p + p0) ? , 


(2.3) 

and the convection tensor 



y _ ^ 

C = p u u , 


(2.4) 

the vector conservation- law form of the momentum equation can be written 

as 

^ (pu) + div T = 0 , 


(2.5) 

where the flow tensor T is defined as 



T = C - P . 


(2.6) 

In Cartesian coordinates x 3> equations 

(2.1) and (2.5) are 


written as 



Cp V + xq (pu 2 5+ ^q Cpu 3^ 

= o , • 

(2.7) 

3 r , 3T 11 „ 3T 12 + 3T 13 

9t Cp V + 8xj + 3 x 2 + 8x 3 

= o , 

(2.8a) 

3 . , 8T 21 3T 22 3T 23 

3t (pu 2 ) + + * xq~ 

= o , 

(2.8b) 

3 , „ 3T 31 3T 32 3T 33 

sr Cp V * 3xj &q~ 

= 0 . 

(2.8c) 


As written, the continuity equation and the three components of the momentum 
equation all consist of sums of derivatives - called a strong conservation-law 
form [2]. The equivalence in structure of these four equations is 

only valid in Cartesian coordinates. In curvilinear, orthogonal coordinates. 
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with scale factors h 2> hj, and unit base vectors ? 2 , e 3 , equations 
(2.1) and (2,5) become 


JT (h 2 h 3 pu l ) * H7 (h l h 3 pU 2^ * (h l h 2 pU 3 D = ° * 


(2.9) 


it + '^z j l 1 T nV + + 4 (hlh2 j =i T ^" ]) = ° 


( 2 . 10 ) 


The decomposition of (2.10) into three scalar equations would introduce un- 
differentiated terms (weak conservation- law form. [2]) due to the terms 

9 e . / 9 x . . Actually, a natural scalar decomposition of (2.10) in the strong 
J ^ 

form is possible [2]. Since this paper is concerned with the formulation 
of the conservation equations as a single vector equation in a higher- 
dimensional space, the vector form of the momentum conservation law (2.10) 
will be retained. The continuity equation (2.9), which is a scalar 
equation, is thus clearly distinct in character from the momentum equation. 


2. 2 Four-dimensional formulation 

The time derivative in (2.10) has an appearance similar to that of the 
spatial derivatives. This suggests that the momentum equation can be put 
completely in divergence form if one considers time as a fourth coordinate 
in a four-dimensional space. Letting x 4 = t, consider the equation 

div 4 T=0, (2.1D 

where T is now a four-dimensional flow tensor. In terms of the components 
of ■?, equation (2.11) can be written as 


_3 

9x 


-(h 2 h 3 h 4 T al t a ) ♦ ^0'l h 3 h 4j 1 + 4; Ch l h 2 h 4.J 1 T a 3 ? a> 


* H7 (h l^ h ^l T a4®a ) - °- (2 ' 12) 
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Greek indexes will indicate quantities defined over the four -dimensional space, 


while Latin indexes will be reserved for quantities defined over the physical, 
three-dimensional space. The independence of time and physical space leads to 
the conditions 


"e . ? = S'e i /3x 4 = 3h./3x 4 = 0 (i «'l to 3). (2-13) 

It follows that £ 4 is constant, so that x 4 is a straight coordinate, and that 
h 4 is at most a function of x 4 » One can thus uncouple the equation for e 4 m 
(2,12). The equation for the remaining components is 


^i-( h 2 h 3 _2 T^?.) + T. 2 e.) 

i 3 - 1 ^ J ~ 


1 3 


+ ^ h i h 2 j f 1 W + W 

Comparing (2.14) with (2.10), one obtains agreement if 


= 0, 


(2.14) 


T j4 " P u j h 4 G- 1 t0 3) i 


C2.15) 


where h 4 is now constant . 

Since the three-dimensional flow tensor has symmetric components T^, 
it is reasonable to assume that the components of the four-dimensional flow 
tensor, , are also symmetric. Therefore let 


T.. = ph.u. 

43 4 3 

The e 4 component of (2.12) thus becomes 

3 3 ~ ' 3 


(j = 1 to 3). 


(2.16) 


a « h.hjl 3T 44 

^C h 2 h 3 p V + + + -pr Sxf " °- (2 ' 173 

One quickly recognizes (2.17) as the continuity equation (2.9), provided 


that 

3T 44 /3x 4 * 0 * (2.18) 

Conditions (2.15), (2.16) and (2.18) can be shown to follow logically 
from the assumption that (2.4) defines a four-dimensional convection tensor 
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?, with the fourth(time) component of the four-dimensional velocity t? satisfying 

u 4 = h 4 . (2,19) 

Equations (2.15), (2.16) and (2.18) can then be replaced by the relation 

T , = T, = C , (a = 1 to 4). (2.20) 

a4 4a a4 v * 

If (2.6) is assumed to be valid in the four-dimensional space, it follows 
that 

P a4 * P 4a = 0 (a = 1 t0 4 >. f2 ' 21) 

Relation (2.19) is made plausible by considering the total derivative as a 
four-dimensional convective derivative, i.e., 

D _ 3 5 U j 3 4 ^ 3 (2.22) 

Dt 3x. . . HT 3x. " . f-i 

4 j = 1 j j a=l "a a 

Alternatively, extending the definition of velocity as the total derivative 
of the position vector to four dimensions, one obtains 

Dx 3x 3 u. 3x 

u a * h a BT - h « W: + ca = 1 t0 4) ' C2 - 23) 

4 J=1 J 0 

Equation (2.19) follows from (2.23) by setting a = 4. 

The analysis, so far, does not provide a natural constant velocity to 
assign to u 4> It can be obtained from a different interpretation of the 
velocity. In terms of the conservation laws, mass and momentum are considered 
fundamental. The velocity can then be defined as the specific momentum, or 
momentum per unit mass. The components of the four-dimensional conservation 
equation in the three spatial directions give the conservation of momentum, 
while the component in the fourth (time) direction has been shown to yield 
the continuity equation expressing the conservation of mass. It is thus 
reasonable to define the fourth (time) component of velocity as the specific 
mass, or unity. 

In summary, it is natural to introduce time as the fourth, straight 
coordinate in a four- dimensional Euclidean space, with the scale factor 
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and velocity component in the fourth direction satisfying 

U 4 = h 4 = i. (2.24) 

The conservation of mass and momentum are then embodied in the single four- 


dimensional vector equation 


div 4 T*- 0, 


with the four- dimensional flow tensor T given by 

-a* -*► 

T * c - P 

where the four-dimensional convection tensor components are 

C aB = pu a u 6 Ca > 8 = 1 to V, 


generali 

zed stress 

tensor 

S u -p-p<fi 

S 12 

. S 13 

S 21 

s 22 -p-p<(> 

S 23 

' S 31 

S 32 

S 33 -p-i 

0 

0 

0 


( 2 . 11 ) 

(2.25) 

(2.26) 

(2.27) 


The in (2.27) are components of the three-dimensional extra stress 
tensor S. Thus, for a constant density flow, a single vector equation 
expresses the conservation of mass and momentum. It is because the time 
coordinate is always straight, and orthogonal to the spatial coordinates, 
that the component in the time direction can be uncoupled to produce a scalar 
equation for the conservation of mass. Since the spatial coordinates can 
be curved, the momentum equation in conservation- law (or divergence) form 
must in general remain a vector equation. 

2 . 3 Steady-State Flow 

An important special case is that of steady-state flow, in which x^ is 
an ignorable coordinate, i.e., 

3/3x 4 = 0 . (2.28) 
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There is a fundamental difference between this case and that in which a spatial 
coordinate is ignorable. The simplest example of the latter is plane flow, 
in which the velocity component in the ignorable direction is zero. The 
component of the momentum equation in the ignorable direction then vanishes 
identically. The number of equations is reduced by one, so that the remain- 
ing equations can be embodied in a single three-dimensional (two space and 
one time) vector equation for a three-dimensional flow tensor. 

In steady-state flow, the velocity component u^ in the ignorable 
direction does not vanish, but remains equal to one. Therefore, the 
component of the four- dimensional conservation equation in the ignorable 
direction does not vanish, but actually remains unchanged. Only in the 
spatial components are the time derivative terms absent. While the number 
of independent variables has been reduced by one, the number of equations 
remains the same. One must still use a four-dimensional formulation to 
describe the conservation laws. While the coordinate is now a fictitious, 
ignorable coordinate from the steady-state viewpoint, the direction of the 
fourth coordinate is very real, since it gives rise to the continuity 
equation. It is appropriate therefore to refer to the coordinate as the 
time coordinate, but to think of the direction of the fourth coordinate as 
the mass direction, since it is connected with the conservation of mass. 

3. CONSERVATION EQUATIONS FOR VARIABLE-DENSITY FLOW 

3. 1 General Formulation 

The continuity equation for variable-density flow is 

+ div (pit) =0* (3.1) 

— 

In the momentum equation, the external body force per unit Volume , pf , can 
no longer be put in divergence form, even if is conservative. It there- 
fore cannot be included as part of a generalized stress tensor, but must 
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appear as an external source terra. The momentum equation is thus written 
as 

(pu) + div t - pf = 0, (3. 2) 

where (2.4) and (2.6) are still valid, but P is defined as 

"? = ? - p?, (3.3) 

The body force ? is left unspecified. In addition to a gravitational force, 
it can also include the fictitious force due to the employment of a non- 
inertial reference frame. In orthogonal curvilinear coordinates, (3.1) and 
(3.2) are written as 


7F Ch l h 2 h 3 P} + H^" Ch 2 h 3 P V + + HJ (h l h 2 pU 3 3 = °> 


(3.4) 


pU jV + dr (h 2V, + T j2«j) 

j-i 1 J * 4“ 3 ■** 

+ ^t h l h 2 T ,S«j) - h l h 2 h 3 P A f j*j ■ °> 

3 J -1 J -1 


(3.5) 


The presence of the variable density requires an additional conservation 
equation to complete the system. The next two sections will treat two dif- 
ferent forms of this additional equation. The results obtained in this 
section will be independent of the specific form of the additional equation. 
The discussion of steady-state flow in ^2.3 showed that a scalar conservation 
equation (in this case the continuity equation) can be associated with a 
coordinate direction, even though the coordinate itself is ignorable. This 
suggests that the additional equation be associated with the direction of a 
new, ignorable coordinate x^.. It thus must have the property that 


5/9x 5 = 0 (3.6) 

for all physical variables. It is easy to show that e,. must be constant, 
and orthogonal to all the other base vectors. Without any loss of generality, 



one can let 


(3.7) 
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The expression for u<- depends on the identity of the additional equation. 
In the five-dimensional conservation equation 

div $ = 0, 

the equations for and can be uncoupled, leaving as the remaining 
equation 

3 3 3 

0 . r* m ^ 1_ r. m s S s* r. m ** ^ 


57]' (h 2 h 3.f 1 T jlV + T j2 e p + ^• (h l h 2 j J 1 T j 3 e j 


(3.8) 


= 0. 


(3.9) 


+ ^ (h l h 2 h 3.f 1 W + ^ h l h 2 h 3.f 1 T jS e ? 

In comparing (3.9) with (3.5), one notes that the source term can be formally 
generated if one lets T. ^ =-pf.x^. Thus the coordinate can be considered 
a source coordinate. This role for the coordinate seems appropriate, since 
the same phenomenon (variable density) requiring an additional conservation 
equation, and an associated new coordinate direction, also necessitates a 
body force term to be treated as an external source of momentum. Now T. r is 
undefined within an arbitrary additive fuction of x^ , x^, x^, x^. Generali- 
zing the concept of the convection tensor to five dimensions, one would 
expect the term pu^u^. In order to satisfy the additional conservation 
equation, an additional arbitrary function gj (x^, X 2 , x^, x^) will be 
necessary. Therefore, let 


T jS = P U j U 5 - pf j X 5 + gj 


() = 1 to 3)j 


(3.10) 


where the are as yet unspecified functions of x^ , X£, x^, x^. 
pu.u. = pu., as in 4 2.2, equation (3.9) is in agreement with the 
equation (3.5). 

If one generalizes (2.20) to five dimensions , i.e., lets 

T 4a = pu 4 u a = pu a (a “ 1 t0 5) > 
then the e^ component of (3.8) becomes 


With T.. = 
momentum 


(3.11) 


+ + iW + = °' ( 3 - 12 ^ 


This is seen to be identical to the variable-density continuity equation (3.4). 
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Assuming that the five-dimensional flow tensor is symmetric, and letting 

3 

g« E g t (3.13) 

j-1 3 M 

the 0^ component of (3.8) takes the form 

div (pu^u) - x 5 div (pf) + div g + (pu^) + ^ — = (3.14) 

Equation (3.14) must be identified with the additional conservation equation. 

The next two sections will consider two such equations of technical interest. 

3 . 2 Incompressible Stratified Flow 

One case of a variable-density flow important in oceanography is that 
of an incompressible, but stratified fluid. The incompressibility condition 
is stated as 

div vl = 0^ (3.15) 

which is equivalent to the conservation of volume. Following the interpretation 
of a velocity component as the specific value of the conserved quantity in the 
corresponding component of the conservation equation, it follows that should 
be defined as the specific volume, i.e., 

u^ = 1/p, (3.16) 

Equation (3.14) can then be reduced to (3.15) if g is set equal to zero, and 
if one defines T^ as 

T 55 = pu 5 + div M x s /2 • (3.17) 

The five-dimensional flow tensor? appearing in (3.8) is thus given by 

T = C - P ■ (3.18) 

where the convection and generalized stress tensors are defined by 

(a,3 = 1 to 5) ; 


C D = p u u D 

aB a 3 


(3.19) 
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and 



S 11 ’ P S 12 S 13 

S 21 S 22 ~ P S 23 

S 31 S 32 S 33 ‘ P 

0 0 0 

pf l X 5 pf 2 X 5 pf 3 X 5 


0 

0 

0 

0 

0 


pf l X 5 

pf 2 X 5 

pf 3 X 5 


div (p?) xjT/2 


(3.20) 


together with equations (2.24), (3.7) and (3,16). 

3. 3 Compressible Equi librium Flow 

The more common case of a variable-density flow is that of a compres- 
sible gas. If the flow is in equilibrium, the energy equation is the only 
additional conservation equation. Let e be the specific internal energy, so 
that the total energy per unit volume is 

e = p [e + j (u 2 + u 2 + u|)] , (3,21) 

Define a heat-flux vector cf to include both conduction and thermal radiation. 

For completeness, let 6 be the net rate of addition of energy per unit volume 
from a non-thermal source, such as a nuclear reaction. The conservation of 
energy can then be expressed as 

+ div (eu + pu - u*S + q) - pu*f - e - 0 , (3,22) 

The velocity component u^ defined by this conservation law is 


u 5 = e/p s e + j (Uj + u 2 + up # 


(3.23) 


Equation (3.14) reduces to (3.22) if one defines 


and 


g = pu - u-S + q j 

T 55 = P u 5 + div X V 2 ” + 6) x 5 . 


(3.24) 

(3.25) 


The five-dimensional formulation of compressible, equilibrium flow is thus 


defined by (3.8), (3.18), (3.19), (2.24), (3.7), and (3.23), with given 


by 
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a$ 


11 


- P 


21 

'31 

0 

i 

15 


12 

; 22 “ P 

S 32 

0 

P 25 


J 1 3 
S 23 

S 33 ' ? 
0 

P 35 


JiVki - pu i - + pf i x s 

k=l u k S k2 - PU 2 ' *2 + P Vs 

k=l u k S k3 ‘ pU 3 ' ^3 + pf 3 JC 5 

0 

( gil U k f k + V V div (P?) 4 /2 


This completes the formulation of the conservation equations for 
variable- density flow. The results will be generalized in ^ 4. 

4. GENERALIZATIONS 


(3.26) 


4.1 Non- Equilibrium Flow 

If the flow of a compressible gas is out of equilibrium, in general one 
must introduce n independent non- equilibrium variables to specify the state of 
the gas. Let such a variable be designated by a , where represents the mass 
fraction of species v, or the specific internal energy (per mass of mixture) of 
the vth internal mode of some species. Equations (3.1), (3.2), and (3.22) are 
still valid, provided that i? is the global velocity of the mixture, and the 
definitions of e and if are appropriately modified to take into account the 
internal molecular structures and interdiffusion of the various species. 

The rate equation for non-equilibrium variable can be written in 
conservation form as 

+ div [pa v (u + u^)] - p * v = 0 (v = 6 to n + 5), (4.1) 

-*d 

Here u^ is the diffusion velocity for a variable which pertains to a 
species in a multi-component mixture, and is the net rate of production 
of due to internal processes. One can then introduce a fictitious, 
ignorable coordinate x^ for each non- equilibrium variable a^, and again 


h^ = 1 {v *= 6 to n + 5), 


define 


(4.2) 
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The velocity component u defined by (4.1) is 


u = a 
v v 


(v = 6 to n + 5), 


The "e component of the (n + 5) -dimensional conservation law 


div r- T = 0 
n+5 


(4.3) 


(4.4) 


is 


nx [ ^7 (h 2 h 3 T vl } + 4r (h l h 3 T V2 :) + ^7 (h l h 2 T v3 :)] 

12 3 1 ° 

3T . 9T _ n+5 9T 

+ + i 'ill =0 (v = 6 to n + 5), 

° x 4 dx 5 3=6 3 

In orthogonal coordinates, equation (4.1) is written as 

9 (pa ,) i % a 3 d 

5t + %^ [h 2 h 3 pa v (u l + U vl\ + ^ h l h 3 pa v (u 2 + u v2 } 

+ ^T^ h l h 2 pa v (u 3 + u v3^* “ p % = 0 (V = 6 to n + 5). 

3 


(4.5) 


(4.6) 


If one again defines 


T v 3 “ p u v u 3 ' P v3 


(v = 6 to n + 5, 3 * 1 to n + 5)^ (4.7) 


then the following relations follow from a comparison of (4,5) and (4.6): 


p . = P . =-n a u . (i = 1 to 3) 

vj jv H V VJ u 


P V5 = P 5v " p Vs 


V ■ p ev ' 0 


(3 ji 1,2, 3, 5) 


(v = 6 to n + 5) , 


(4.8) 


4. 2 Magnetohydrodynamic Flow 

If the fluid is electrically conducting, it is subject to external 
forces and energy transfer due to an electromagnetic field. For non- 
relativistic fluid motion, displacement currents and charge accumulation 
can be neglected, and Maxwell’s equations are 
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-j* 

div b =» 0 f 


(4.9) 


9b 


^ + curl e = 0 ; 


(4.10) 


curl b - Hj = 0^ 


(4.11) 


where b, e, j are the magnetic induction, electric field, and current density, 

respectively, and H is the constant magnetic permeability of a vacuum. 

- ^ 

The external force per unit volume, j X*>, can be rewritten, with the 
aid of (4.9) and (4.11), as 


j X b = ~ curl b x b = ± div (bb - J—L. i) . 


(4.12) 


The energy transferred to the fluid per unit volume from the electromagnetic 
field, j*e, becomes 


“* 1 — * 1 — V — * IS o 

j*e = jj-curl b*e = - - div (eXD ^ (|b| 2 ), ( 4 . 13 ) 

The effect of the electromagnetic field is to add to the five-dimensional 
generalized stress tensor P^g the electromagnetic generalized stress tensor pj^ 
defined by 



b l ' I z b k 
1 k=l K 

b l b 2 

b l b 3 

0 

b 2 e 3 ' 

b 3 e 2 

b 2 b l 

2 3^2 

b 2 ‘ 7 z b k 

Z Z k-1 K 

b 2 b 3 

0 

b 3 e l - 

b l e 3 

b 3 b l 

‘ b 3 b 2 

2 1^2 
b 3 • \\ b l 

k=l 

0 

b l e 2 - 

b 2 e l 

0 

0 

0 

0 

1 3 
'J 1 

k 2: 

b k 

b 2 e 3 ‘ b 3®2 

b 3 ? l - b l e 3 

b l e 2 - b 2 e l 

3 

1 y h 2 
’ 2 1 b k 
k=l 

k=l 

0 



Note that the ten, independent non-zero components of P m Q can be identified 

Ctp 

with the components of the four-dimensional, symmetric Maxwell stress-energy 
tensor defined in the relativistic formulation of electrodynamics. 
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4.3 General Coordinates 

For many important applications, one must use coordinate systems more 
general than orthogonal coordinates. Let the physical space be described by 
general, curvilinear coordinates xi, X 2 , This defines a system of covariant 

base vectors ^ = Br^/Sx 1 , where is the position vector. The covariant 
components of the metric tensor are then given by g ^ The formulation 

of the conservation equations in general coordinates will be illustrated by 
considering compressible equilibrium flow. One therefore introduces a straight, 
fourth coordinate = t, and a fictitious, ignorable, straight, source 

r 

coordinate x . The corresponding constant covariant base vectors, which are of 
unit length, are ^ and The covariant components of the five-dimensional 

metric tensor are thus given by 



*11 

g 12 

g 13 

0 


g 21 

g 22 

*23 

0 


g 31 

*32 

g 33 

0 


0 

0 

0 

1 


0 

0 

0 

0 


Let the determinant of.g^ be designated by g = Ig^l 
The five-dimensional conservation law 

div^ T = 0 

can be written in component form as 



(4.15) 


(3.8) 


A t-E ** ■ «, 


(4.16) 


where summation over repeated indexes is assumed, Greek indexes being summed 


from one to five, while Latin indexes are summed from one to three. The 

a8 

contravariant components of the flow tensor, T , can be expressed as 
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ja$ = c a3 _ p a3 ( a ,3 = l to 5)^ 


(4.17) 


where the convection tensor is defined by 

C a6 = pu a u 6 (a, 6 = 1 to 5) (4.18) 

with 

u = 1 ; (4.19) 

and _ . . 

u = e/p = e + j uV, (4.20) 

Here u 1 (i = 1 to 3) are the contravariant components of the real, physical 

velocity. The contravariant components of the five-dimensional generalized 

stress tensor take the form 


,ot3 


s 11 - 

11 

Pg 

s 12 - 

12 

Pg 

s 13 - 

13 

Pg 

0 

j c il 
g ij uS 

1 

- pu - 

1 

q 

+ 

J 5 

pf X 

o 21 

21 

c 22 

22 

c 23 

23 

0 

g..u j S i2 

2 

2 


r 2 5 

S 

Pg 

s 

Pg 

5 

Pg 

- pu - 

q 

+ 

pf X 

S 31 - 

31 

Pg 

S 32 - 

32 

Pg 

s 33 - 

33 

Pg 

0 

g..uV 3 

-T 

3 

“ pu - 

3 

q 

+ 

_3 5 
pf x 


0 

.15 


0 

.25 


P 35 0 


(pg..u*f^ + e)x 3 - — (^pf *)—,.— 

2 


) 


(4.21) 

where g 1 * 1 , S 1 ^, q 1 , f 1 are the contravariant components of the metric tensor, the 
extra stress tensor, the heat-flux vector, and the external body force per unit 
mass, respectively. 

The extension to general coordinates of other forms of the conservation 
equation is carried out in an analogous manner. 


5. DISCUSSION 

The multi-dimensional conservation equations (2.11), (3.8), or (4.4) do 
not constitute a complete set of equations. These equations involve two classes 
of dependent variables. One class consists of the primary variables. For a 
compressible gas, these are u^, p, and e. For non- equilibrium flow they also 
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include a^. For an incompressible fluid, e must be replaced by P* The other 
class consists of the secondary variables.' The latter include S ; ^ (or 

40 , e, and p (for a compressible gas). For non- equilibrium flow, and 
must also be introduced, while for magnetohydrodynamic flow the field component 
and b^ appear as secondary variables. All of these secondary variables must 
be related to the primary variables in order to obtain a complete system of 
equations. 

The relations between the two classes of variables may involve the intro- 
duction of a third class of intermediate variables. Expressions for the heat- 
flux vector components normally involve the temperature. The radiant heat- 
flux is most generally calculated in terms of the specific intensity. The 
electromagnetic field components are calculated from Maxwell's equations, which 
involve the current density vector components j^. In addition, particular forms 
of constitutive equations for the extra stress tensor, heat-flux vector, dif- 
fusion velocity, and current density require the introduction of various trans- 
port coefficients such as viscosity coefficients, diffusion coefficients, thermal 
conductivities, electrical conductivities, absorption coefficients, etc. 

The particular forms that constitutive equations take, as well as 
expressions for equations of state, transport coefficients, and the internal 
production terms (5^, are all based on assumed models for a given fluid. These 
models lie outside the scope of the basic formulation of the conservation 
equations for the fluid motion. One can, however, categorize these models as 
to their effect on the completeness of the set of conservation equations. For 
a simple fluid model the constitutive equations give the secondary variables 
as explicit functions of the primary and intermediate variables, and their 
derivatives. The multi-dimensional conservation equations then constitute a 
complete set of differential equations (supplemented by algebraic relations 
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for the intermediate variables) in the absence of electromagnetic fields, and 
in the presence of thermal radiation treated within the Rosseland diffusion or 
the emission-dominated approximation, (In the latter case, the divergence of 
the radiant heat -flux can be treated as an external energy source term, 6.) 

For a complex fluid model, the expressions for the secondary variables 
could involve integrals of the primary variables, or require the introduction 
of new variables which are given by auxiliary equations. Thermal radiation 
generally requires the solution of the radiative transfer equation. It is 
coupled to the flow equations via the absorption, emission, and scattering 
coefficients. The presence of electromagnetic fields requires the solution of 
Maxwell's equations, which are coupled to the flow equations via a generalized 
Ohm’s law. In all of these cases, the conservation equations form only part 
of the complete system of equations describing the physical situation. 

As an example of the application of the multi-dimensional formulation, 
consider the compressible equilibrium flow equations written in general 
coordinates ($4.3). The most general coordinate transformation in- 

volving the three physical coordinates and time is 


ci 1 1 2 3 4 V c 

x a = x a (x ,x ,x J ,x 4 ) x* = x 5 (a’ =1 to 4) f (5.1) 

4 « 

where x is the new time like coordinate. In order to achieve a simple 
form of scalar decomposition, base vectors and tensor components will be 
defined with respect to the unprimed coordinate system. The transformed 
Eq. (4.16) thus becomes 


a. 1 

3 trT 3x a 


t 6v = 0 . 


(5.2) 


Since g^ and ^ are constant, one obtains immediately the transformed 


continuity and energy equations 
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9x 


a 1 




T 3x 


a' 


3x 


v 


T® V ) = 0 


ce = 4,5) 


(5.3) 


(5.4) 


1 ' 2 T 3' 

The remaining base vectors can be functions of x , x , x an 
x 4 ’. In any finite-difference numerical procedure, all unknown variables 
are defined at discrete points fixed in x , x ,x ,x space. Let the 
subscript P’ denote a quantity evaluated at one of these points. A 
subscript a ! appearing to the left of a quantity will denote variation with 
x at only, the other coordinates being held fixed. at their values at point 
P*. Introduce unit base vectors by 

g (i) " ^ y/g ii ' ; " 

In any finite-difference algorithm applied to (5.2) at point ?' , the unit 
base vector g (j) at a point other than P« , appearing in the x a derivative 

term, can be expressed as 

■f - V M (5.5) 

a' g (j) a' Y (j) g (k)P’ . 

The expansion coefficients a .Y^j are known functions of x“ , determined by 
the original coordinate system and the transformation (5.1). The scalar 
decomposition of the transformed momentum equation in strong conservation 
form becomes 

3x“ 3x w_ a ' y W 

If one specialized (5.3) and (5.6) to an ideal gas (S = q = f = 6 = 0), and 
introduces physical contravariant velocity components by 

^ 5 /g77 u 1 (unsummed), 


L_ ^gTTT je ) = 0 (k .1 to 3) , 

a’ 18 ZT a ,T 0) *i) 


(5.6) 


u 


(5.7) 


one obtains the following set of strong conservation equations: 


a a* 4 ’ 3x 4 ' uCU 


:)] + “rr [*'5 7 ‘ PC^ft - + , j , — 

Sx 1 


■i » 


3x* U 


(j) 


:)] - o. 


n 


(5.8) 
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3x 


3 i f—r .,00 r3x^ (j) . 3x^ ^pu^u^ ; jm,,^ 

r W 4 ^(j) h*- pu + — y-— ' ■ + P^T g )]} 


(5.9) 


3x 


^mm 


3 r r~r „,(k) r^x" (j ) 3x ^pu^u*" ^ } j^Kil n 

+ — W i,Y iri hST pu + — C ~ "^ ' + g )]} K 0 


3x 


3x 


(k = 1 to 3), 


4’ 


(j) 


3x 


ri*f t 


3x J 


’JJ 


(5.10) 


♦ -tt' f/g 1 " e + = o. 

Sx 1 3x J Jg 77 . 

S J1 

Further details and specializations of these results may be found in 
Ref. [2]. 


Another application of the higher-dimensional formulation is the 
calculation of fluid motion in a non-inertial reference frame fixed with 
respect to a moving body. In the conventional approach, external body force 
terms appear in the momentum equation, and an external work term is 
introduced in the energy equation.. One can show that the change in 
reference frame can be treated as a coordinate transformation in a higher- 
dimensional space, with the external source terms arising from the 
curvature of the new time coordinate. Using a method analogous to the 
one just outlined, one can then obtain a scalar decomposition of the 
equations in strong conservation- law form. The details will appear in 
a future publication. 
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